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ABSTRACT 


The object of this work is to formulate an analytical 
model of the heat transfer phenomena in deep submergence 
submarine shaft seals, It is expected that this will prove 
useful to design engineers in determining materials and 
arrangement specifications for these seals. 

Interfacial temperature data derived from a comprehen- 
sive experimental investigation is compared with thermal 
results computed from a lumped model of the shaft seal. 

This latter analytical data is obtained by first defining 
the model, then calculating the heat source component in- 
puts, and finally utilizing an IBM 7094 computer to aid in 
solution of the resulting mathematical entanglements. 

It was found that the most significant factor in the 
overall heat transfer scheme is the quantity of bypass cool- 
ing water'flow. No bypass flow results in an extremely sharp 
temperature rise with pressure; whereas only a slight amount 
of flow will reduce the interfacial temperature rise to ac- 
ceptable levels. It was also found that with a combination 
of high inlet pressure and low shaft speed, the largest tem- 
perature rises were experienced. 

It can be concluded that any design change that would 
limit or eliminate the source of cooling water must be care- 
fully evaluated as to the thermal implications. Further, in 
any investigation of this nature, due consideration ought to 
be given to the low-speed/high-pressure limiting condition. 
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INTRODUCTION 


Limitation*tof the amount of sea water that penetrates 
the hull of a submarine along the propulsion shafting pre- 
sents a variety of problems, most of which nave been or 


×0 ٦ 


are being examine However, one area of interest, 
in which no experimental work has been Gone, and little ana- 
lytical effort has been expended, is in the thermal limi- 


2 
(2,3) mne type of shaft seal 


tations of submarine shaft seals. 
presently being utilized in submarine design is of a radial 
face, balanced pressure design with a small amount of leakage 
flow. (see Figure I) A requirement exists to determine the 
quantity of heat generated at the seal interface and the 
efficiency with which this heat can be dissipated. 

Up until quite recently, no heat generation or heat 
dissipation problems were encountered in the design of these 
seals. This was due to the fact that submarines were designed 
with an auxiliary salt water piping system in which the pres- 
sure was maintained slightly above sea pressure, It was only 
a secondary purpose of this system to provide positive flow 
to the shaft seal. This provided a slight amount of leakage 
through the seal and a large amount of bypass flow overboard. 


This latter bypass flow served to effectively remove generated 


heat from the seal. Two conditions have arisen which now make 
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the thermal limitations of these seals a problem that must be 
faced and quantítatively determined. First, the auxiliary 
salt water system may be completely eliminated in future 


(6) 


submarine designs. This would effectively remove a 
major portion of the necessary heat sink of the shaft seals. 
Further, as submarine designs and operatiug profiles trend 
more and more towards deeper submergence, the quantity of 
heat generated at the seal face grows yguite rapidly. 

1t is the object of this thesis to first present a 
description of an experimental investigation into the ther- 
mal conditions existing in deep submergence submarine shaft 
seals. Then an analytical model of shaft seal is devel- 
oped and utilized to formulate a method of predicting maxi- 
mum interfacial temperatures, In the analysis, some of the 
variables are determined solely from theory, while the de= 
termination of others is based in part on theory and partly 
on the results of the experimental data. The experimental 
results are then presented and compared with the analytical 
results. Extrapolations are made to include conditions be- 
yond the physical capabilities of the test rig. 

It is expected that the analytical procedure determíned 
in this paper will be utilized by shaft seal design engineers 
in material and arrangement specifications, and that a great 
deal of costly and time consuming experimental work will 


thereby be avoided. 





II 


EXPERIMENTAL PROCEDURE 


A 5 5/8" shaft diameter test rig for submarine shaft 
seals was utilized for gathering experimental data, Maxi- 
mum testing conditions were limited to a pressure of 2000 
psi and a rotational speed of 375 rpm. 

A bearium seal ring was selected for instrumentation. 
Bearium (commercial trade name) is a leaded bronze compo- 
sition that has been found quite suitable for use as a 
seal ring in submarine shaft seals. (see Figure II) Five 
1/64" holes were drilled at equal (0.044") spacings across 
the 0.263" of the seal ring face. These holes were coun- 
ter-sunk and fitted with 1/8" copper tubing for further 
connection to pressure sensing devices. (see Figure III) 
Holes=were drilled for nineteen (19) chrome alumen thermo- 
couples which were then inserted in the seal ring and seated 
with epoxy. (see Figure IV) 

The seal ring face was then lapped and profiled with 
an X-Y plotter at three separate circumferential positions. 
The face ring (tungsten carbide) and the seal ring were 
inserted in the housing and assembly of the test apparatus 
was completed. (see Figures V and VI) The area ratio (face 
area of the sealing surface to balance area behind face ring 


housing) was 0.523. A sealol step seal, for which leakage 
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FIGURE Մ 


SEAL RING INSTALLED IN FRESSURE HOUSING 
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FIGURE VI 


TEST SEAL AND CONTROL SEAL INSTALLED IN HOUSING 
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and torque data was already available from previous testing, 
was inserted in the back half of the test rig. Instrumenta- 
tion was set up to measure the five (5) pressure taps, nine- 
teen (19) thermocouple temperatures, torque, and rpm. (see 
Figure VII) Filtered river water (brackish) was utilized 

as the working fluid. Data was taken with the following 


fixed variables: 


p N Q, 
200 psi 50 rpm 0.0 gpm 
200 " oO met 0.1 " 
500 ۳ 29 ն 1,٢۴ 기 

1000 " 375 ٥ 
1500 EN 
2000 " 


Besides the instrumented measurements, inlet pressure, inlet 
water temperature, seal leakage flow, and seal leakage flow 
temperature were measured and recorded. 

Upon completion of all data taking, the seal ring face 
was again profiled at the same locations. The results of 


these profiles (before and after) are shown in Figure VIII. 
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FIGURE VII 
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III 


ANALYTICAL PROCEDURE 


Before determining the quantity of heat generated at 
the seal interface and the various parameters which are in- 
volved in that operation, it was first deemed appropriate 
to try to set up an analytical model of the seal apparatus. 
Various techniques were considered and it was finally decided 
that a lumped model based on thermal resistances and a gross 
heat source was best suited for this investigation. Thermal 
resistances were carefully determined based on the most 
accurate predictions of heat coefficients available and 
on judicious estimation of the areas and linear distances 


(8) 


of the relatively complicated seal geometry. The con- 
vection quantities were determined by the experimentally 
fixed bypass flow and measured seal leakage flow. (see 
Appendix A) The final model as utilized throughout the 
analysis is displayed in Figure IX. 

The overall heat input at the seal face is a combina- 
tion of a variety of heat source components which vary in 
magnitude with the physical parameters. Basically, there 
are three calculable contributors to the total heat genera- 
س88‎ ٦ These are the heat generated by the dissipation 


of pressure energy, the heat generated by viscous shear, 


and the heat assocíated with boundary lubrication. Thus 
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the total heat generated at the sealing interface can be 


calculated from the following equation: 


4 = Anress * Ivise * bound 


OT, 


2 
Q هم‎ MAU fAC, UADp 


00۶ ۷ً 7 

In the above equation, only the quantity of flow through 
the sealimg faces "Qe" is an experimentally measured quan- 
tity. Pressure drop "ap" and rotational speed "U" are fixed 
variables, while face area "A" is fixed by geometric consider- 
ations and "sg" and "J" are constant terms. Therefore, those 
variables remaining to be determined are the coefficient 
estimating the net average contact stress "Cp"? the seal 
face clearance "e", and the coefficient of friction "f", 

The coefficient estimating the net average contact stress 
E can be calculated by utilizing the results obtained 
from the five (5) pressure taps installed in the seal face. 
The pressure ranges recorded at the five taps spaced at 
radially equal distances apart are plotted in Figures X, 
XI, XII, XIII, XIV, and XV. The range of pressures measured 
for various running conditions of "N" and ues is plotted 


against a non-dimensional constant "K", where "K" is the ratio 
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FIGURE XII 
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FIGURE XV 
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of the difference in radii of the outer seal face and the 
pressure tap to the radial distance between pressure taps. 
Now, since "Co is a function of ,"م"‎ and the pressure 


profiles are known, the seal unbalance may be pictured as: 


Dp —> 


LS d ose [E 
| 
| 
| 





_ Closin ressure - openin ressure (3) 
p inlet pressure 


For example, at "p = 2000 psi", (see Figure X): 


0 . 0.52 - (p/6) — 6+0.05+ 0.02+0.01) (4) 


0.523 = 0.249 = 0.274 

The balance pressure thus determined and the seal face leakage 

flow are shown plotted against inlet pressure in Figure XVI. 
The seal face clearance "c" for viscous radial flow 

with "e!" and "w" constant can be determined by summing the 

pressure and shear forces across the fluid film and applying 


Newtonts Law of víscous flow. In simple form: 


2 
A 4d u 
ar B (5) 


Since the velocity "u" is equal to zero at "y=0" and "y=c", 
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FIGURE XVI 
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a yo = cy) (6)‏ و 
Finally, since‏ 


Գբ 


Շ 
այ Udy; (7) 
Q = TER . (8) 


In equation (8), "c" was determined by using a weighing 
function on "w", based on the percent of seal face actually 
experiencing the opening pressure, This weighing function 
was determined from the pressure profiles of Figures X thru 
XV. For "p = 1000 psi", "2/3 w" was assumed; for "p = 200, 
300, and 500 psi", "5/6 w" was assumed; and for "p = 1500 
and 2000 psi", "1/3 w" was assumed. Seal face clearance 
is shown graphically in Figure XVII. 

The determination of a reasonable coefficient of fric- 


tion "ft was finally resolved in the following manner: 


t t (in-1b) (9) 


Ն րթ» 20.2 Cop (psi) 


Torque data was available from two sources. Original ex- 
periment produced some data, but difficulty with the mechan- 


ics of the rubber coupling at which the torque-meter was 
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attached made the results somewhat dubious. Data from sim- 
ilar tests by K. S, Saste111 ^! provided the author with 
another source of plotted torque evidence. Both of these 
Sources were weighed and the combined most accurate predic- 
tion of torque was determined, The coefficients of friction 
thus calculated are shown plotted in Figure XVIII against 


a duty parameter "G" where 


«նս -9 N (rpm 
G = = 1.224 x 10 «ՅԵՐ (10) 


ք 


With all the variables now accounted for, the total 
heat input at the seal face source was calculated for all 
combinations of the fixed parameters: 

p 200 = 2000 psi 


N 150 - 375 rpm 
ԳԵՏ 0.0 = 1.0 gpm 


From Figures XVI, XVII, and XVITI, the variables were 
extrapolated to include pressures of 3000, 4000, and 5000 psi 
in order to display the effect of very deep submergence on 
the heat input and interfacial temperatures. Total heat 
input, determined as described on the previous pages, is 
shown graphically as a function of inlet pressure and shaft 
speed in Figure XIX. 

With the heat input as source, and Figure IX as model, 
an IBM 7094 computer was utilized to determine the nodal 


temperatures of the resulting 5x5 matrices. (see Appendix B) 
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IV 


RESULTS 


The following thirteen (13) graphs make up the body 
of thermal experimental and analytical results of this 
investigation, In Figure XX is shown the actual seal face 
temperatures (determined analytically) versus the source 
heat input. Figures XXI and XXII show actual seal face 
temperature (determined analytically) plotted against inlet 
pressure for bypass flows of 0.0 and 0.1 gpm respectively. 
Figures XXIII through XXXII depict both an analytical and 
an experimental non-dimensional temperature quantity "T", 


where 


CTI) 


This is plotted against nz the bypass flow rate, for 
varying inlet pressure from 200 to 5000 psi. Note that 
Figure XXVIII is a duplicate presentation of Figure XXIX 
with the exception of the scale utilized. Note also that 
all data in Figures XXX, XXXI, and XXXII is from the analyt- 


ical extrapolation. 
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FIGURE XX 


SEAL FACE TEMPERATURE VS. HEAT INPUT 
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FIGURE XXI 


SEAL FACE TEMPERATURE VS, PRESSURE 
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FIGURE XXII 
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V 


DISCUSSION OF BESULTS 


Ihe first three figures of Chapter IV present, in cap- 
sule form, the major points learned from this investigation. 
The most obvious and also most important result is best seen 
in Figure XX on page 30. With no bypass cooling water flow, 
the seal interfacial temperature rises sharply with heat in- 
put, so swiftly in fact that a heat input of only about 
1600 BTU/HR is required for vaporization of the film.  How- 
ever, it can be seen that only a very slight quantity of 
bypass cooling water flow (0.1 gpm) will decrease the tem- 
perature rise to a safe level for even the highest predicted 
value of heat input. Further, if the quantity of bypass 
cooling water flow is increased to 1.0 gpm, the seal inter- 
facial temperature rise is almost negligible for even large 
increases in heat input. 

Figures XXI and XXII on pages 33 and 34 respectively 
display the same phenomena when examined together. Taken 
separately, some further discussion presents itself. In 
Figure XXI, it can be seen that for zero bypass cooling water 
flow, the limiting seal inlet pressure is between 3400 psi 
and 3700 psi, depending on the rotational speed of the ap- 
paratus. A very interesting phenomenon also discernable in 


this Figure is the sharper rise in temperature of the low 
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speed analytical data at high values of input pressure, 
This can be attributed first to the results of Figure XIX 
on page 28, which in turn is really a function of the three 
heat input variable terms. Further examination of the ana- 
lytical data shows that the key lies in the heat associated 
with boundary lubrication, the third term in the heat equa- 
tion (Equation (2) on page 15). For example, the low value 
of duty parameter "G=1.34x107 108 at "p=5000psi" and "N=150rpm" 
results in a relatively high value of the coefficient of 
friction "f=0.018", as can be seen on the extrapolated por- 
tion of Figure XVIII on page 27. 

The same phenomenon is also discernable (to a lesser 
extent) in Figure XXII on page 32. Here, with a small amount 
of bypass cooling water flow, the temperature rises are not 
very sharp, and are also not significantly dependent on RFM. 

At very low pressures, the spread of results between 
the experimental and the analytical data was a maximum. (see 
Figures XXIII, XXIV, and XXV) However, at higher values of 
the seal inlet pressure (p = 1000 psi and above) the experi- 
mental and analytical results show a marked similarity in 
trends and values, It is extremely interesting to note that 
the trend for the combination of low RPM and high pressure 
to result in a greater heat input and hence higher inter- 
facial temperature begins at a considerably lower pressure 


in the experimental results then it does in the analytical data. 
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The fact that this trend did exist both experímentally and 
analytically is quite significant. 

The curves of Figure XXVIII on page 38 display in true 
focus the sharp rise in the temperature associated with de- 
creased bypass cooling water flow. It would have been de- 
sirable to display all of the temperature increases in this 
form, but this would have resulted in the curves being so 


bunched as to be unreadable; hence the logarithmic displays. 
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VI 


CONCLUSIONS AND RECOMMENDATIONS 


Deep submergence submarine shaft seals of a similar 
design to those presently being used require a minimum 
amount of bypass cooling water flow. Any design change 
that would limit or eliminate the source of cooling water 
must be carefully evaluated as to the thermal implications 
involved, it is recommended that the method of evaluation 
should follow that presented in the Analytical Procedure 
and in the Appendices of this paper. 

Further, in any investigation of this type, due consid- 
eration should be given the low-speed/hish=-pressure limiting 
condition. This warning is especially emphasized since one 
would intuitively expect that a high-speed/high-pressure 
limiting condition would be all that required examination. 

In future experimental investigations with submarine 
shaft seals, it is recommended that a more satisfactory method 
of gathering accurate torque data be devised and utilized. 
Only then can a great deal of confidence be given the computed 
coefficients of friction and hence heat source inputs. 

Finally, computer programs similar to that described in 
Appendix B are available in most large computation centers, 


and ought to be utilized in investigations of this nature. 
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APPENDIX A 


DETERMINATION OF THERMAL BESISTANCES 


The thermal resistances shown diagrammatically in Fig- 


ure IX on page 14 were calculated from the following basic 


equations: 
q = h AAT 2 (12) 
1 ոճ 
_XK د‎ ° Te UI 
d + © CT Hum ce (14) 
9 Q e 


The following constant coefficients were established: 


Bearium 20.4 
Monel 1550 
Tung. Carb, JOTO 
Epoxy 0.6933 


Cooling Water ل‫ د0‎ 





Monel-air 
Interface 
Bearium-Monel 
Interface 

Monel-Water 
Interface 
Bearium-Water 

Interface 
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The value of "Q," was set either at 0.0 (closed cir- 
cuit) or at 0.1 or 1.0 gpm. Values of "Qe" were determined 
individually for the range of pressures utilized. (see Fig- 
ure XVI on page 23) 

The following cross-sectional areas of heat transfer 
surface and lengths of heat transfer paths were estimated 
by examination of the particular geometry of the seal under 
study. A similar analysis would be required for any modi- 
fication in seal geometry. However, for a simple scaling 
modification, the areas and linear distances could be simi- 


larly scaled. 


Ու ոշ) | ը (ո) | 


1.5 


Bearium x 31.0 — | 
Bearium-Monel Interface | (TO | 
Monel (seal ring housing) | 7 | 192 
Monel-Air Interface | 1100 

(seal ring housing) | | 
Bearium-Water Interface | 34.0 | 
Sealing Water at Face | 271 | 1.4 
Sealing Water at Leakoff 1 0 | ee | 
Seal Water-Air Interface | CD | 
Tungsten Carbide տ. Oe 
Epoxy | 34.0 | 0.005 
Monel (face ring housing) | 22,0 120 
Monel (shaft) | 150.0 10 
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Monel-Air Interface 
(shaft) 

Monel-Water Interface 
(face ring housing) 

Bypass Water at Face 
Ring Housing 

Water-Monel Interface 
(seal housing) 

Monel (seal housing) 


Monel-Air Interface 
(seal housing) 
Bypass Water at 
Bypass Channel 
Bypass Water-Air Interface 
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APPENDIX B 


COMPUTATIONAL PROCEDURE 


The computer program utilized in this thesis can be 


found in the Department of Naval Architecture and Marine 


Engineering Computer Library at the Massachusetts Institute 


of Technology under the coding of NA 021.0. The computer 


listing follows: 


500 
400 


501 
502 
503 


دم دا دک د OO‏ 


00 


16 


LO 
1 
12 
5 


FORTRAN 

DIMENSION A(11,11),B(11,11),X(11), E( 11), IL( 20) 
READ 400, JOB 

FORMAT (I5) 

IF (JOB) 501, 502, 503 

CALL EXITM 

CALL EXIT 

IF (JOB-2)100,200,300 

PART 1 SIMULTANEOUS EQUATION SOLN BY XSIMEQF 


F FACTOR BY WHICH XSIMEQF MULT. VALUE OF DET. A(I,J) 
JJ-NUMBER OF COLUMNS IN MATRIX B 
M=1 SOLN GOOD,BRIGHT BOY,A(1,J) REPLACED BY X(1I,J) 
M=2 OVERFLOW, BROKE THE BANK 
M=3 MATRIX A IS SINGULAR, GO TO JAIL 


READ 15 

PRINT 15 

READ 10,N,((A(I,J),I=1,N),J=1,N) 
PRINT 11,N 

DO 16 I=1,N 

PRINT 19,(A(I,J),J=1,N) 

CONTINUE 

PRINT 24 

READ 2,(B(K,1),K=1,N) 

PRINT 2,(B(K,1),K=1,N) 

FORMAT (5F13.4) 

PRINT 17 

PRINT 18 

PRINT 24 

FORMAT (15/(5F13.4)) 

FORMAT (4H N= ,13) 

FORMAT (16H GOOD BOY DET = ,1PE14.6/(op5F14.5)) 
FORMAT (16H BROKE THE BANK) 
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14 


15 


18 


19 
24 


25 
26 


20 
21 


22 
23 


FORMAT 

LECT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


(54H GO TO JAIL, DO NOT PASS GO, DO NOT COL 
200 DOLLARS) 

( 54H ) 
(1H2) 

(8H ANSWERS) 

(1HO,10F11.4) 

(180) 

(F10.4) 

(8H SCALE= ,F10.4) 


77606257٦ 
PRINT 26,F 


PRINT 
آ37‎ 


շկ 


MEXSIMEQF(11;N,JJ,A,B,F,E,) 


GO TO 
PRINT 
GO TO 
PRINT 
GO TO 
PRINT 
GOSTO 


(2021,22) IN 
12,F,(A(I,1),1=1,N) 


The program, when supplied with the proper input control 


cards, 


will solve simultaneous algebraic equations up to a 


maximum of 10 equations and 10 unknowns: 


a. When the solution is successful = BRIGHT BOY. 
Ե. When overflow occurs - BROKE THE BANK, 
c. When the matrix is singular = GO TO JAIL. 


d., Whe 


n matrix to be inverted is not square - 
(caustic comment). 


A complete write-up of this FAP program can be obtained from 


the MIT Consultants (cc-174=3). 


INPUT 


Card 1 


Format (1 


This is the control card. Place a 1 in column 5. 


Card 2 Format (54H ) 


Flace 


a 1 in column 1 of the data card, then type 


Po 





any comment that you wish to have printed. 

Card Format (I 
The number of unknowns in the set of simultaneous 
equations is punched in columns / and 5; the units 
digit appearing in column 5. 

Group 4 Format (5F13.4) Cards = دود‎ 
The value of the coefficients are punched five to 
a card. The coefficients are punched column by col- 
umn (i.e., first column - all rows, second colume - 
all rows, etc.). 

Group 5 Format (5F13.4) Cards = 5 
The value of the algebraic equations are punched 
flve per card. There is no need to fill a11 5 
fields on the card. 

Card 6 _ Format (F10.4) 
This value serves as a scale factor by which XSIMEQF 
multiplies the value of the determinent of the mat- 
rix, If the value is one, a 1.0 is punched anywhere 
in the above field. The scale factor is printed out 
in the output data as SCALE = „ This scaling option 
applies to those who are manipulating large numbers. 

Sample input and output data follow. The comment card 

(1st card) indicates in order: inlet pressure (psi), bypass 


flow (gpm), and rotational speed (rpm). 
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200 T10 
N= 5 

576.40 

~ 35.00 

210.00 

-18,40 

~ 513.00 


50520 


ANSWERS 

SCALE= 

GOOD BOY 
087.206 


2/35 


-35.00 
1436.72 
-1400.00 
Or 
O. 


148.00 


1:00 


10.00 
-1400.00 
1420.00 
0: 
o 


860,00 


DET = 3.619504E 13 


-18.40 
O. 
O; 
15973109 
-1500.00 


1609.00 


519200 
0. 
0. 
=1500.00 
2529.00 


4400.00 


CX GHO ap > GP OD ap QP A GD աթ AD P cp AP Ap سه‎ Q Q Q Q GO Q GP COND O Q G AGED = G Q QD Q ae CUP GD a os E) GES SF GP s Q Q Q Qs ab 


11000 ۶) 


N= 5 
1045 
-35.00 

-6.05 
-18.40 
= 51.40 


205.70 


ANSWERS 
SCALE= 
GOOD BOY 


92.84 


87.48 87.47 
150 
=35.00 -6.05 
1436.72 -1400.00 
-1400.00 1412.05 
0. O. 
O. 0. 
148.00 516.00 


1.00 


DET = 6,280575E 11 


DIE 


005ھ 


-18.40 
OF 
OF 
او٢‎ 
-1500.00 


1609.00 


-51.40 
0. 
er 

-1500.00 
1605.20 


4629.99 





3000 
N= 5 
54,24 
35.00 
-0.72 
-18.40 
=0.12 


106. 70 


ANSWERS 
SCALE= 
GOOD BOY 


۱۵ 0 


0 250 


- 0 
1436.72 
-1400.00 
0% 
0. 


148.00 


. 00 


163 25 


-0.72 
-1400.00 
1401.41 
0. 
0. 


59.30 


0 10 


70700۲ 


-18.40 
0 . 
0. 
1537.10 
=1 500,00 


1609.00 


-0.12 
03 
0 


-1500.00 
1502.72 


224,00 


In all, ninty-one (91) data points (solution of 91 


separate 5x5 matrices) were obtained utilizing this pro- 


cedure, 


The five temperatures given in the output of 


this program conform respectively to the five lettered 


(A,B,C,D,E) nodes in Figure IX on page 14. 
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